The heavy quark potential in pNRQCD 
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The heavy quarkonium static potential is discussed within the framework of potential NRQCD. Some quantita- 
tive statements are made in the kinematical situation ran 3> Aqcd at the level of accuracy of the next-to-leading 
order in the multipole expansion. 



1. INTRODUCTION 

In a recent series of papers ||]|] a detailed study 
of a suitable Effective Field Theory for heavy 
quark bound states, called potential NRQCD 
(pNRQCD) ||, has been started. Since several is- 
sues have been treated in that context, we address 
the reader to the quoted literature for a complete 
overview of the achieved results. Here we only 
mention the 1-loop matching and the static 
energies of the hybrids S presented in these pro- 
ceedings. While in the following, due to its con- 
siderable importance, we recollect and summarize 
our understanding of the heavy quarkonium static 
potential. 

Let us define, first, what we mean with heavy 
quarkonium potential. Being heavy quarkonium 
a non-relativistic bound system, it is character- 
ized by at least three energy scales: the mass or 
hard scale, m, the momentum or soft scale, mv, 
corresponding to the inverse of the bound-state 
size and the energy or ultrasoft (US) scale, mv 2 
(v is the heavy quark velocity). As a consequence, 
when US degrees of freedom are neglected, heavy 
quarkonium can be described as a bound state 4> 
governed by a non-relativistic Schrodinger equa- 
tion of the type 



^ + U(r,p,S 1 ,S 2 ,m) )<t> = E<j>. 
m 



U) 



We will call V the heavy quarkonium potential, 
which is in general a function of the quark dis- 
tance r, momentum p, spin Si, S2, and mass m. 



Another relevant quantity for heavy quarko- 
nium physics is expected to be the energy between 
static sources. For heavy quarkonium in a singlet 
state this can be defined as 



E s {r) = lim -ln(W n ) 



(2) 
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where Wn is the static Wilson loop of size r x T 
and the symbol ( ) means the average over the 
gauge fields. Often in the literature E s has been 
implicitly identified with the static limit of the 
Schrodinger potential, lim V . While it is rea- 

m — >oo 

sonable to expect this identification to hold to 
some extent, there are no general grounds for it 
to be true in general. Indeed, already long time 
ago H several doubts have been rised on the in- 
frared consistency of that identification at least 
in perturbation theory. It is the goal of this con- 
tribution to make some quantitative statements 
on the difference E s — lim V. This will be done 

m — >oo 

in the last section. In the next section we per- 
form, as an intermediate step, the matching in 
the singlet sector of the pNRQCD Lagrangian at 
the next-to-leading order in the multipole expan- 



2. pNRQCD 

Another scale is relevant in QCD, the scale 
where nonperturbative effects start to become im- 
portant. We will call this scale Aqcd and we 
will assume that mv 3> Aqcd- For sufficiently 
heavy quarkonium v -C 1 and therefore the en- 
ergy scales of the system are widely separated. 
This allows to systematically integrate out these 



scales by matching QCD with simpler but equiv- 
alent Effective Field Theories. The integration of 
the hard scale (~ to) gives rise to the effective the- 
ory known as non-relativistic QCD (NRQCD) 
whereas the integration of the soft scale (~ mv) 
gives rise to what we call potential NRQCD. Be- 
ing m and mv well above Aqqd both matchings 
can be done perturbatively. 

By definition pNRQCD is the effective the- 
ory where only degrees of freedom below the soft 
scale remain dynamical. The surviving fields are 
quark-antiquark states (with US energy) and glu- 
ons with energy and momentum below mv. It 
is convenient to decompose the quark-antiquark 
states into singlets and octets under colour trans- 
formation. The relative coordinate r, whose typ- 
ical size is the inverse of the soft scale, is explicit 
and can be considered as small with respect to 
the remaining dynamical lengths in the system. 
Hence the gluon fields can be systematically ex- 
panded in r (multipole expansion) . Therefore the 
pNRQCD Lagrangian is constructed not only or- 
der by order in I/to, but also order by order in 
r. As a typical feature of an effective theory, all 
the non-analytic behaviour in r is encoded in the 
matching coefficients. 

The most general pNRQCD Lagrangian den- 
sity that can be constructed with these fields 
and that is compatible with the symmetries of 
NRQCD is given at order l/m° (but we write also 
explicitly the kinetic energy in the centre-of-mass 
frame) and at the leading order in the multipole 
expansion by: 



£ p nrqcd = Tr|s f ^id - ^~ 



-Ot [ iD - — - V {r) 

TO 



V s (r) 



O 



+gV A (r)Tr {O f r • E S + S f r • E 0} 
+.9-^^Tr{O t r -EO + OtOr-E} , (3) 

where S = S(r, R, t) and O = 0(r,R, t) are 
the singlet and octet wave functions respectively 
and R is the centre-of-mass coordinate of the 
quark-antiquark system. All the gauge fields in 
Eq. (0) are evaluated in R and t. In particular 



E = E(R, t) and iD O = id O - g[A (R, t),0]. 

We call V s and V the singlet and octet static 
matching potentials respectively. By looking at 
the equations of motion of the Lagrangian (||) it 
is clear that, as far as higher order terms in the 
multipole expansion (terms of order r or smaller 
in (||)) do not give potential-type contributions, 
V s and V coincide with the static singlet and 
octet potential to be used in the heavy quarko- 
nium Schrodinger equation. This happens when 
the US scale mv 2 is the next relevant scale of the 
system (i.e. Aqcd < mv 2 ). While, in the situ- 
ation mv 3> Aqcd 3> m v 2 one expects to have 
nonperturbative corrections to the static poten- 
tial coming from higher order terms in the multi- 
pole expansion. Both situations will be discussed 
in the next section. 

Here we sketch the singlet matching at order 
1 /to and at the next-to- leading order in the mul- 
tipole expansion. We refer the reader to [0J|] for a 
complete and detailed discussion. The matching 
is in general done by comparing 2-fermion Green 
functions (plus external gluons at a scale below 
mv) in NRQCD and pNRQCD, order by order 
in 1/to and order by order in the multipole ex- 
pansion. In order to get the singlet potential, we 
choose the following Green function in NRQCD: 



I = S 3 ( Xl -yi)5 3 (x 2 -y 2 )(IU D ), 



(4) 



where Wn is the rectangular Wilson loop with 
edges x x = (T/2,r/2), x 2 = (T/2,-r/2), y x = 
(-T/2,r/2) and y 2 = (-T/2, -r/2). In pN- 
RQCD we obtain at the next-to-leading order in 
the multipole expansion 



Z s (r)5\^ 



yi)<5 3 (x 2 - y2 )e-^« 

F/2 dt f dt'g-itt-OO'.-v.) 

-T/2 J-T/2 



(r.E a (t)«KM'):" j r-EV)> 



(5) 



where 4> J is a Schwinger (straight-line) string 
in the adjoint representation and fields with only 
temporal argument are evaluated in the centre- 
of-mass coordinate. Comparing Eqs. (f|) and (||), 
one gets at the next-to-leading order in the multi- 
pole expansion the singlet wave-function normal- 



ization Z s and the singlet static potential V s . Va 
and V must have been previously obtained from 
the matching of suitable operators, but for the 
present purposes we only need the tree-level val- 
ues: Va = 1 and V a = (Ca/2 — Cp)ct s /r. Let us 
concentrate here on the matching potential V s . 
By substituting the chromoelectric field correla- 
tor in Eq. (||) with its perturbative expression we 
obtain at the next-to-leading order in the multi- 
pole expansion and at order In a s 



V s (r)=E s (r)\ 



2-loop+NNLL 



r it 12 2r[i 



(0) 



where E s has been defined in Eq. (|^). We note 
that V s and E s would coincide in QED and that 
therefore the effect we are studying here is a gen- 
uine QCD feature. The 2-loop contribution to E s 
has been calculated in ||. The NNLL contribu- 
tions arise from the diagrams studied first in Q 
and shown below. An explicit calculations gives 



C F C\at Xn C A a s 



12irr 2 
-0(1/T) 



Inserting this and the 2-loop contribution in Eq. 
(f§ we get 



V s (r) 



a v (r, /j,) 



a v (r, = a s (r) < 1 + (ai + 2^ E Po) 



+ 



16 7T 2 



a s (r) 

47T 



7£ (4a 1 /? + 2/3 1 ) + ( — +4 7 |)/3 2 



+a 2 



C\al{r) 

12 7T 



In r/i > , 



(7) 



where f3 n are the coefficients of the beta function 
(a s is in the MS scheme), and a\ and a 2 are given 
m§. 

We conclude this section noticing that the octet 
matching potential V D can be calculated in the 



same way as done for the singlet. In this case the 
relevant NRQCD Green function could be chosen 
to be 5 3 ( Xl - yi)5 3 (x 2 - y 2 ){T a W n T a ), where 
the colour matrices are inserted in the endpoint 
Schwinger strings. Even if this Green function is 
gauge dependent, as discussed in g], the match- 
ing should guarantee a gauge invariant definition 
of V . A 2-loop calculation is still not available, 
but work is in progress 

3. THE STATIC SINGLET POTENTIAL 

In the previous section we have established 
the connection between E s , the energy of static 
sources and V s the singlet static matching poten- 
tial of pNRQCD. Here we discuss in two different 
kincmatical situations the connection of V s with 
the static limit of the heavy quarkonium potential 
defined through the Schrodinger equation. 

A) Aqcd < m» 2 . This situation is expected 
to hold for toponium and for the bottomonium 
(charmonium?) ground state. As already men- 
tioned, in this situation V s , as given by Eq. (0), 
is the heavy quarkonium static potential in the 
sense given in the introduction. The explicit /i de- 
pendence of it originates from the fact that the US 
degrees of freedom (which have the same scale of 
the kinetic energy and therefore do not belong to 
the potential) have been explicitly subtracted out 
form the static Wilson loop. This fact is not sur- 
prising if we understand the heavy quarkonium 
potential as a matching coefficient of pNRQCD. 
As a consequence even in a purely perturbative 
regime the static heavy quarkonium potential (as 
well as ay) turns out to be an infrared sensi- 
tive quantity. In this situation nonperturbative 
effects are only of non-potential nature (see for 
instance the Leutwyler-Voloshin type corrections 
in the situation Aqcd *C mv 2 |1C]]). Finally we 
mention the quite obvious fact that, when calcu- 
lating any physical observable, the \i dependence 
in (0) must cancel against /Lt-dependent contribu- 
tions coming from the US gluons (see for instance 

0)- 

B) mv ^> Aqcd 3> mv 2 . Since in this situation 
there is a physical scale (Aqcd) above the US 
scale, a potential can be properly defined only 
once this scale has been integrated out. At the 



next-to-leading order in the multipole expansion 
we get 



V(r) 



-C F 



a v (r, fx) 



r 



g 2 i, y ' x 



dte-iUVo-v.) 
x(E a (O0(t,O)^E b (O))(/i), (8) 



e -«(V -V.) = 1 _ it( y o _ Vs) _ _ Fs) 2 + 



Therefore, the heavy quarkonium static poten- 
tial V is given in this situation by the sum of 
the purely perturbative piece calculated in Eq. 
(0) and a new term carrying also nonperturba- 
tive contributions (contained into non-local gluon 
field correlators). This last one can be organized 
as a series of power of r n . We stress that due to 
the condition mv 3> Aqcd this expansion makes 
sense only in the short-range. Typically the non- 
perturbative piece of Eq. (^) absorbs the /i depen- 
dence of ay (see 0] for an example) so that the 
resulting potential V is now scale independent. 

The infrared sensitivity of the static potential 
can also be expressed in terms of renormalons 
(see for instance Q). Rephrasing them in the 
Effective Field Theory language of pNRQCD we 
can say that the singlet matching potential V s , as 
defined in Eq. ([?]), suffers from IR renormalons 
ambiguities with the following structure 



V s {r)\ 



IR ren 



Co 



C 2 r 2 + 



(9) 



The constant Co ~ Aqcd is known to be 
cancelled by the IR pole mass renormalon 
(2m polc | mrcn = -Co, ||). While Eq. (|) pro- 
vides us with the explicit expression for the op- 
erator which absorbs the C% ~ Aq CD ambigu- 
ity More precisely the order r 2 term on the 
right-hand side of Eq. (j|) suffers from UV and IR 
renormalons. The UV renormalon ambiguity of it 
(which can be calculated simply by substituting 
the chromoelectric field correlator with its per- 
turbative expression and summing up the leading 
log of all the bubble diagrams) cancels exactly the 



second term in the expansion (^): 



(10) 



x di(E-(t)0(t ) O)^E 6 (O))(/ i ) 
Jo 



= -C 2 r 2 



UV ron 
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